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In this work we present a detailed explanation of the construction of an appropriate equation 
of state (EoS) for nuclear astrophysics. We use a relativistic model in order to obtain an EoS for 
neutrally charged matter that extends from very low to high densities, from zero temperature to 
100 MeV with proton fractions ranging from (no protons) to 0.6 (asymmetric matter with proton 
excess). For the achievement of complete convergence, the Sommerfeld approximation is used at 
low temperatures and the Boltzman distribution for relativistic particles is used in the calculation of 
the electron properties at very low densities. Photons are also incorporated as blackbody radiation. 
An extension of this EoS is also presented with the inclusion of strangeness by taking into account 
g ■ the E" hyperon only. Strangeness fractions range from 0.02 to 0.3. 

o 

^ ■ I. INTRODUCTION 

sa 

1 The observed pulsars, also commonly known as neutron stars, are believed to be the remnants of type II supernova 
^*xplosions. These type II supernovae appear at the end of the evolution of very massive stars. The core of these stars 
V^collapses to a density around several times nuclear saturation density. A rebound takes place and drives a Shockwave 
O^hich cxpells most of the original mass of the star. The simulation of supernova explosions and the conditions for 
them to take place have been subject of investigation for the last 30 years. Depending on certain thermodynamical 
T T^onditions present in the equations of state (EoS), the supernova explosion simulation is successful or not [1]. The 
i^»EoS built for nuclear astrophysics purposes depends on a series of thermodynamic properties which are obtained 
■^-obr certain temperatures, densities and matter composition. Hence, an efficient EoS which is reasonably accurate is 
v s Qnandatory for a supernova explosion simulation to be successful. 

C^) 'In order to obtain an equation of state (EoS) for low and high density matter suitable to astropysical applications, 
Qbe relativistic non linear Walecka model (NLWM) HQ is used. For matter to be neutral, electrons are also included. 
ttTFor sufficiently high densities the formation of hyperons is energetically favored. Normally, the appearance of the 
^^trange baryons softens the EoS. Our formalism is described next with the inclusion of the whole baryonic octet for 
nfFp sake of completeness but, in a first step towards a complete desciption of a supernova explosion, only protons and 
Oieutrons are considered. Next we incorporate strangeness but restrict ourselves to the inclusion of £~. Convergence 
^problems are well known to exist at low temperatures and below certain densities. Appropriate approximations are 
5-4tlen utilized. Blackbody radiation is considered and, whenever convenient, electrons and positrons are described 
separately. Future prospects for obtaining more sophisticated EoS are discussed. 



II. HADRONIC MATTER EQUATION OF STATE 



A common extension of the NLWM considers the inclusion of the whole baryonic octet (n, p, A, E + , S°, S~, 
E ) in the place of the nucleonic sector. In this work we present the complete formalism, but numerical calculations 
were performed with nucleons only. 

The lagrangian density of the NLWM reads: 

£ ~\~ ^mesons ~T" ^leptons 3 (1) 



where 



C B = [l P - 9vbW - g pB t ■ W) - {M B - g sB <j))] yj t 



with extending over the chosen baryons B, 

9sB = XsB 9s, 9vB = x vB 9v, 9pB = x pB 9p 

and x s b, x v b and x p b are equal to 1 for the nucleons and acquire different values in different parametrizations for 
the other baryons, 

Cmesons = \{d ^ <j> - m^ 2 ) - 1 K 3 - ^Xcj> A - + 
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-iB^.B^ + im^-lA (2) 

where f2 M „ = d jL V v — d„V^, B MI/ ~ d^h u — <9„b p — <? p (b M x b^) and t is the isospin operator. 

In the above lagrangian, neither pions nor kaons are included because they vanish in the mean field approximation 
which is used in the present work and we do not consider the possible contribution of pion and kaon condensates. The 
leptonic sector is included as a free fermi gas which does not interact with the hadrons. Its lagrangian density reads: 



Cieptons = y^V'i (il^ -n%i)ipi. (3) 
i 

In the present work only electrons (and positrons are considered). In the mean field approximation (MFA) (see 
[3,0], for instance), the meson equations of motion read: 

00 = -7r^4>l - :A^o +y2^2 X ^B PsB, (4) 

2m% 6mj z ^ mj 



Vo=y~] —^XvB PB, (5) 
= £ ^ X PB *3fl PB, (6) 



with 

PB 



2 J r A(fB + ~fB-), P = £PB, (7) 



1 ( A/f* 
PsB = -r / P 2 dp , = UB+ + /*-), 



7T^ 



P 2 + Mg 



with A/^ = Mb — g S B <j>, B± stands respectively for baryons and anti-baryons, t^s is the third component of the 
baryon isospin, E*(p) = \/p 2 + M* 2 and 

/ S ± = l/{l + exp[(£T(p)T^)/T]}, (8) 

where the effective chemical potential is 

vb = Pb - 9vbV - g pB t 3 B b . (9) 

Within the MFA, the meson fields are taken as classical fields whilst the baryon fields remain quantum |2j. On the 
other hand, the Dirac equation, which is the equation of motion for the baryons is not solved directly but instead 
used in the calculation of the densities appearing in the meson equations of motion. The system has then to be solved 
self-consistently. 

At T = 0, the distribution functions for baryons are replaced by step functions. In this case equation becomes 
simply pb = kp B /3ir 2 . The baryonic energy density in the mean field approximation reads: 



f *8 + + ±4 + ^ + ^ + ^bi (io) 



and the related pressure becomes 



3tt 2 ^ 7 Vp 2 + M* 2 
ml 2 /c^g A$ to 2 , 2 ™L 2 

"2™ 6 24~ ~ + — b °- (U) 
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The entropy of the baryons are taken as 

l3 



Sb = - 2 T.j 4^lfB+ l °9(fB + ) + (l-fB + )log(l-f B +) 

+ f B -log(f B -) + (1 - f B -)log(l - /b-)] (12) 

and hence the free energy reads 

T B = £b - TS B - (13) 

Notice again that the above expressions were obtained for finite temperature, but they can be easily modified for 
T = 0. Whenever T=0, no anti-particles are present. 

At very low T (—0.4 < log(T) < —0.1 MeV) there are well known convergence problems due to the distribution 
functions and in this case we use the Sommerfeld approximation for the baryons [6|. The effective chemical potentials, 
in particular, read 

Vi = e Fi - — T ±- '-, i=p,n. (14) 

o KFiC-Fi 



For the net electron density we have 



/'■ = -> /7CT(/.- "A*). ( 15 ) 



where the distribution functions for the particles (e ) and antiparticles (e + ) are given by 

/ oT =l/(l + exp[(e T A*e)/T]) ) (16) 
with fi e as the chemical potential. In order to ensure charge neutrality, electron and proton densities have to be equal, 



i.e. 



Pe = P P - (17) 

Next we always distinguish between electrons (e~) and positrons (e + ) and when both particles and antiparticles are 
considered we refer to the related quantity with the index e. At T = 0, the distribution functions for the leptons are 
also replaced by step functions and no positrons are left. In this case equation l|15p becomes simply p e = kp e /3ir 2 . 
The thermodynamic quantities read 

£ e = 2 / ^P^T^ e (fe- + /«+), (18) 

E e - =2 / ^VP 2 +m] f e -, £ e+ =2 [ ^LPy/pi+ml /e+j (19) 
J (zir) J \2it) 

P e = -L f P * dp (/ e -+/ e+ ), (20) 

3?F •/ VP + He 

P '' /y ' /e", P e+ =Pe-P e - (21) 



3tt 2 7 ^p 2 + m 2 e 



_ E e + P e - H e p e _ £ e - + P e - - [i e Pe- Q _ £ e + + -P e + + PePe+ , 
— = , O e — , O e + — — , [ZZ) 



(23) 



4 



The particle fraction is defined as yi = p%/p, where i = p, n, e~ , e + , and p is the total baryonic density. 

At very low densities a Boltzman distribution for relativistic electrons and positrons is necessary Q. The low 
density limit depends on the temperature and is numerically chosen such that eqs. (fl"5l) and (fl7|) are equal within a 
1CP 6 precision. If the difference is larger than this limit, eq. (TIT)) is chosen and the corresponding chemical potentials 
are 

(24) 

with 3 = 2 defined as the spin multiplicity, Moreover, 

Pe = (e^ T -e-^ T )^, (25) 

or analogously, 



Me = Tl0g[~ + ^ + 1], Z = 7T 2 p e /h. (26) 

The energy density and pressure become 

E e = {e^ T -e-^ T )^ (27) 



log 



Pe / 2?r 
9 \Tm ( 



3/2 



Pe = {e ^T_ e -^T ) (h-mll ) ^ 

dir z 

where (3 = 1/T, J = ^W), h = ~f = f^i(y) - f ^, I a = -% = ^(y) - 2$^ + f 
with j/ = m e (3, Ki are modified Bessel functions and K' v (x) — —\ (K u -i(x) + K v+ i{x)). 
The photons are taken into account via blackbody radiation and the main expressions are 



For the hadron phase we have used the GM3 parametrization proposed by Glendenning and Moszkowski [f|, 
corresponding to an effective mass M* — 0.78 M and incompressibility K = 240 MeV at the saturation density 

p a = 0.153 fm- 3 . The coupling constants are = 9.927, {j^j = 4.82, (j£} = 4.79, k = 0.017318gs 3 , 

A = -0.014526 ff s 4 . 

In our codes the inputs are the temperature, proton fraction and baryonic density. The grids for these quantities 
are -0.4 < log(T) < 2 (MeV) with mesh intervals of 0.1, 
< y P < 0.6 with mesh intervals of 0.02, 
3 < log(p) < 15.7 (g/cm 3 ). 

In the output we have y p , y n , y e ~, y e +, y e , 
p p - M (MeV), p n -M (MeV), p e - m e (MeV), -p e + m e (MeV), 
£ B (erg/g), £ e - (erg/g), £ e + (erg/g), £ 1 (erg/g), £ B + £ e + £ 1 (erg/g), 

Sb (^s/baryon), S e - (/cs/baryon), S e + (fcs/baryon), S 1 (fcs/baryon), Sb + S e + S~ ( (fc^/baryon), 
Pb (dyne/cm 2 ), P e - (dyne/cm 2 ), P e + (dyne/cm 2 ), P~ l (dyne/cm 2 ), P B + P e + P-y (dyne/cm 2 ), 
T B +Fe+F~t (MeV/fm 3 ). 



III. THRESHOLD DENSITY FOR MATTER WITH STRANGENESS 



To include strangeness in the EoS, S _ was first chosen because in (3- equilibrium matter at zero temperature and 
with the GM3 parametrization, its onset appears at lower densities than the onset of the least massive hyperon, the 
A. Depending on the parametrization chosen for the NLWM and on the hyperon-meson coupling constants, this trend 
may change at higher temperatures and hence, in a future work £~ and A should be included simultaneously. 

In compact stars, stellar matter is in chemical equilibrium, which means that 



PT.0 — PS" = PA = Pn, MS- = PB.- = Pn + Pe, Ps+ = Pp = Pn - Pe- 
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In an explosive enviroment like the one existing in a supernova, chemical equilibrium is not supposed to be enforced. 
However, we consider that the time during which the supernova explosion occurs is much longer than the characteristic 
time of the weak interaction in such a way that the strangeness fraction is expected to be finite. 

In order to build an EoS containing strangeness and appropriate for a supernova simulation we define for each 
energy density, temperature and proton fraction a threshold density above which a given fraction of strangeness, y s , 
is allowed to exist. We determine the threshold density from the condition of /3-equilibrium for the E - , which is 
imposed through the two independent chemical potentials (fi n and fi e ). In this case, at T = 0, the corresponding 
effective chemical potential and density are 



Vy,- = \/k z FTi + M* s = /i s - - g vS Vo + g P xb , 

and 



Ps 



"■fs 



3?r 



2 ' 



If the condition 



Ps . 

— > Vs, 
P 

is satisfied, the appearence of the strangeness fraction y s in the EoS is allowed. For y s > 0, we define 

P = Pn+ Pp+ PS 

with 

p P = y P ( 1 -ys)p, Pn = 0--yp)0--ys)p, ps = y s p. 

For charge neutrality 

Pe= P P - Ps- (30) 

is required. 

In order to fix the meson-hyperon coupling constants we have used the prescription given in [j| [^], where the 
hyperon coupling constants are constrained by the binding of the A hyperon in nuclear matter, hypernuclear levels 
and neutron star masses (a^s = 0.7 and x v j: = £ps = 0.783) and assumed that the couplings to the E are equal to 
those of the A hyperon. 



IV. RESULTS AND FUTURE PROSPECTS 

A comprehensive test of the thermodynamic accuracy and consistency of our EoS, as described in [Toj |. mainly when 
strangeness is introduced, should be performed. 

A non-homogeneous phase known as pasta phase should be considered at low densities. This non-homogeneous 
configuration made out of spheres, rods, bubbles or other more exotic structures, have been extensively used recently 
These strucutres may change the neutrino opacity in supernova matter and influence neutron star quakes 
and pulsar glitches. We can obtain the pasta phase by building the binodal section, and therefore obtaining the 
chemical potentials and densities of the gas and liquid phase in equilibrium. A very crude approximation would be 
to forget Coulomb interaction and take zero thickness nuclei. We can consider the matter made of liquid dropplets 
in a gas introducing two parameters: the radius of the Wigner-Seitz cell and the radius of the nucleus (equal for 
protons and neutrons). One of the parameters would be fixed imposing a given particle density and the other by 
minimizing the free energy. These results can be improved by including the Coulomb contribution and the surface 
energy by hand. A Thomas-Fermi calculation can then be used to obtain the pasta phase with all fields introduced 
in a consistent way and the surface energy calculated from the derivatives of the fields. 

a-particles can be easily incorporated in the EoS as proposed in [Tl| . Once the pasta phase and the a-particles are 
included, the EoS should then be compared with reference [15j . We are already aware of some important differences. 
The parametrization used in [lij ]. known as TM1 [l6| reproduces ground state properties of stable and unstable nuclei. 
Nevertheless, this parametrization has proven not to be adequate in the description of neutron star matter because it 
breaks down, giving rise to negative baryon effective masses at densities exisiting inside a neutron star (approximately 
6 times the nuclear saturation density) when hyperons are incorporated into the EoS 17]. For this reason, we usually 



G 



choose one of the parametrizations introduced by Glendenning and collaborators [9(, which give a higher nucleon 
effective mass at the nuclear matter saturation density and, for this reason, avoids the problem of the baryon negative 
masses. Moreover, according to [Hj], the EoS with inhomogeneties has a critical temperature T ~ 15 MeV above 
which matter is uniform. This number certainly depends on the choice of the parameters. Based on our recent works, 
we would expect a smaller value for the critical temperature since for nuclear matter with no electrons (no Coulomb 
interaction and surface tension) the critical temperature occurs for symmetric matter just above 15 MeV. The high 
value obtained in [l5| may be due to the way the density distributions are parametrized which give rise to very stiff 
surfaces for the droplets. In [l3| a critical temperature of ~ 5 MeV was obtained for y p = 0.3 matter and ~ 6 MeV 
was obtained for y p — 0.5. One of our recent studies on the dynamical instabilities of npe matter also predicts lower 
critical temperatures, more according to the results of 13]. We do not know if the differences on the EoS due to the 
use of different parametrizations is more important or of the order of the magnitude of the changes included due to 
the explicit inclusion of a non-homogeneous phase. This should be studied. 

In Fig. [1] we compare our results for the free energy obtained at three different temperatures and three different 
proton fractions with the results of fl5l |. One can see that results deviate sligtly at higher densities. In Fig. [2] we 
plot the pressure for the same temperatures and proton fractions as in Fig[TJ Again the results are very similar. In 
Fig. [3] we plot, once more, for the same temperatures and proton fractions, the entropy. The differences are more 
pronounced. While at higher temperatures (50 MeV) the curves are very similar for all proton fraction, for lower 
temperatures the curves are identical only for neutron matter (very low proton fraction). One should notice, however, 
that the trends of the curves are the same. 

Finally, we comment on the definition of the internal energy: it is equal to the nucleon mass for zero density at 
T=0 MeV. For finite temperature this is no longer true because of the presence of nucleons and antinucleons. We 
have defined the internal energy as the energy density per number density in erg/g and in [l5j the internal energy is 
given by the energy density per number density minus the atomic mass unit. As one can see in Fig 21 both results are 
in accordance once the same definition is used. A more clear comparison is done in Fig. [5] where the internal energy 
for homogeneous matter within TM1 is also shown, and compared with the internal energy obtained with the GM3 
parametrization and the EoS of [l5j], also with TM1 but with the non-homogeneous phase included. 

Finally, as a second step in a more refined EoS with strangeness the A hyperons and later the whole octect and 
muons should be included. 



ACKNOWLEDGMENTS 

This work was partially supported by CNPq(Brazil), CAPES(Brazil)/GRICES (Portugal) under project 100/03 
and FEDER/FCT (Portugal) under the projects POCTI/FP/63419/2005 and POCTI/FP/63918/2005. 



[1] J. Cooperstein, H.A. Bethe and G.E. Brown, Nucl. Phys. A 429, 527 (1984). 

[2] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

[3] J. Boguta and A. R. Bodmer, Nucl. Phys. A292, 413 (1977). 

[4] D.P. Menezes and C. Providencia, Phys. Rev. C 68, 035804 (2003). 

[5] S.S. Avancini and D.P. Menezes, Phys. Rev. C 74, 015201 (2006). 

[6] N. Ashcroft, N.D. Mermim, Solid State Physics, Saunders College Publishing, Orlando, 1976 

[7] L. D. Landau and E. M. Lifshitz, Statistical Physics, Pergamon Press, 1959. 

[8] N. K. Glendenning and S. Moszkowski, Phys. Rev. Lett. 67, 2414 (1991). 

[9] N. K. Glendenning, Compact Stars, Springer- Verlag, New- York, 2000. 
[10] F.X. Timmes and D. Arnett, Astrophys. J. Suppl. 125, 277 (1999). 
[11] J.M. Lattimer and F.D. Swesty, Nucl. Phys. 535, 331 (1991). 

[12] D. G. Ravenhall, C. J. Pethick, and J. R. Wilson, Phys. Rev. Lett. 50, 2066 (1983); M. Hashimoto, H. Seki, and M. 

Yamada, Prog. Theor. Phys. 71, 320 (1984). 
[13] G. Watanabe, K. Sato, K. Yasuoka, and T. Ebisuzaki, Phys. Rev. C 69, 055805 (2004); G. Watanabe, T. Maruyama, K. 

Sato, K. Yasuoka, and T. Ebisuzaki, Phys. Rev. Lett. 94, 031101 (2005). 
[14] T. Maruyama, T. Tatsumi, D. N. Voskresensky, T. Tanigawa, and S. Chiba, Phys. Rev. C 72, 015802 (2005). 
[15] H. Shen, H. Toki, K. Oyamatsu, K. Sumiyoshi, Nucl. Phys. A 637, 435 (1998); H. Shen, H. Toki, K. Oyamatsu, K. 

Sumiyoshi, User Notes for Relativistic EoS Tables. 
[16] K. Sumiyoshi, H. Kuwabara, H. Toki, Nucl. Phys. A 581, 725 (1995). 

[17] D.P. Menezes and C. Providencia, Phys. Rev. C 68, 035804 (2003); A.M.S. Santos and D.P. Menezes, Phys. Rev. C 69, 
045803 (2004). 



y p =0.5 



T=50, yp=0.5, ours 
Shen 

T=10, yp=0.5, ours 
Shen 

T=5, yp=0.5, ours 
Shen 



y p =0.3 




T=50, yp=0.3, ours 
Shen 

T-10, yp=0.3, ours 
Shen 

T=5, yp=0.3, ours 
Shen 



y P =0.01 




T=50, yp=0.01, ours 
Shen 

T=10, yp=0.01, ours 
Shen 

T=5, yp=0.01, ours 
Shen 




FIG. 1: Free energy (MeV/fm 3 ) as function of the energy density (g/cm 3 ) for different temperatures (MeV). 
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FIG. 2: Pressure (MeV/fm 3 ) as function of the energy density (g/cm 3 ) for different temperatures (MeV). 
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FIG. 3: Entropy per baryon as function of the energy density (g/cm 3 ) for different temperatures (MeV). 
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FIG. 4: Internal energy of symmetric matter for different temperatures. 




FIG. 5: Internal energy of symmetric matter for different temperatures for homogeneous matter within GM3 and TM1 
parametrizations, and the EoS of [l5l ]. 



